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Abstract. We present the results of CO(1-0) emission mapping with the IRAM interferometer, at ~ 1 " resolution, 
of the z=0.223 ultra-luminous starburst IRAS 11582+3020. This galaxy was selected from an IRAM-30m survey 
of 30 galaxies at moderate redshift (z ~ 0.2-0.6) to explore galaxy evolution and in particular the star formation 
efficiency, in the redshift range filling the gap between local and very high-z objects. The CO emission is kinemat- 
ically resolved, and about 50% of the total emission found in the 27" (97 kpc) single dish beam is not recovered by 
the interferometer. This indicates that some extended emission may be present on large scales (typically 7-15"). 
The FIR-to-CO luminosity ratio follows the trend between local and high-z ultra-luminous starbursts. 
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1. Introduction 

The star formation history (SFH) of the Universe has been 
extensively studied in recent years, through high spatial 
resolution observations of star-forming galaxies as a func- 
tion of look-back time (e.g. Madau et al 1998). The main 
striking feature is a steady increase of the star formation 
rate between z=0 and z=l, by at least a factor 10 or even 
more (Blain et al 1999a). 

Starbursts were more frequent in the recent past, and 
star-forming regions were also more dust enshrouded - the 
evolution of galaxies appears to be much faster in the in- 
frared than in the optical/UV. This strong evolution of in- 
frared bright star-forming galaxies is not really understood 
in models of galaxy evolution. It has been established ob- 
servationally that the number of galaxy mergers, trigger- 
ing starbursts, increases with redshift as a high power- 
law, like (1 + z) 3 (e.g. Le Fevre et al 2000), and this is 
consistent with numerical simulations and semi-analytical 
models (e.g. Balland et al 2003). However to reproduce 
the strong SF evolution, the contribution and efficiency of 
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mergers to the star-formation must also vary considerably 
with redshift, with a peak at z ~ 1, to agree with the 
observations (cf. Blain et al. 1999b, Combes 1999). 

To better understand the physics responsible for this 
evolution, it is paramount to measure the content and dis- 
tribution of the fuel for star formation, the molecular gas, 
as a function of redshift. Although a few tens of ultra- 
luminous infrared galaxies, often amplified by a gravita- 
tional lens, have been mapped in the CO line at very high 
redshift (e.g. Omont et al 2003, Walter et al 2004, Tacconi 
et al 2006), not much is known about star- forming galax- 
ies at moderate distances, in the range 0.2 < z < 0.6 
(only 2 objects have been studied in this range, Solomon 
et al 1997). Therefore, we have undertaken a CO survey 
of 30 IR-luminous galaxies in this redshift range to check 
whether the derived molecular gas content and star forma- 
tion efficiency (SFE or SFR/M(H2)) are evolving at this 
faster rate as well. 

For this study, we first obtained global information 
on the molecular content and on the SFE using sin- 
gle dish data 1 . After a successful detection in our single 
dish survey, we obtained high-resolution Plateau de Bure 
Interferometer (PdBI) observations to map the molecu- 

1 An article presenting the results of our galaxy survey will 
be presented elsewhere, Combes et al. in prep. 



2 



F. Combes et al.: High resolution CO map of a starburst at z=0.223 



lar gas distribution. Resolving the gas distribution and 
kinematics is crucial to attribute the star formation activ- 
ity to a global merger or a less violent process. The first 
high-resolution CO observations of the survey have been 
obtained towards IRAS 11582+3020 which is the object 
of this Letter. 

2. The sample and the source IRAS 11582+3020 

The ULIRG sample of Solomon et al (1997) contains 37 
objects, however only 2 have z > 0.2. To fill the gap be- 
tween low and high-z (z > 2) studies, we have initiated 
the first systematic survey of 0.2 < z < 0.6 sources, se- 
lecting the 30 most luminous galaxies detected at 60 mi- 
cron (IRAS or ISO) with known spectroscopic redshifts. 
IRAS 11582+3020 is one of the 12 well detected galax- 
ies. It is an ultra-luminous galaxy with L(IR)= 5.4 10 12 
L (with an IR flux computed from IR= 1.8 x (13.48 
fi2 + 5.16 f 25 +2.58 feo+fioo) 10~ 14 Wm" 2 ; Sanders & 
Mirabel 1996). Based on the ion lines in its optical spec- 
trum it was classified as a LINER by Kim et al (1998). The 
source is not detected in 2cm radio continuum by the VLA, 
with an upper limit of 0.8 mJy (Nagar et al 2003), but it 
is a 20cm radio-source, with a flux of 3 mJy (Becker et 
al 1995). The FIR-to-radio ratio g=log([F F/i? ./(3.75 10 12 
Hz)]/[f„(1.4 GHz)])=2.7, typical for ULIRGs (Sanders & 
Mirabel 1996). Rupke et al (2005) find evidence for a su- 
perwind outflow in this galaxy of about 15 M yr , while 
its SFR is estimated at 740 M yr _1 from the infrared lu- 
minosity, SFR= ah IR /(5.8 10 9 L o ) (e.g. Kennicutt 1998), 
using a correction factor of a=0.8 for AGN contribution, 
as adopted by Rupke et al (2005). The red image of Kim 
et al. (2002) reveals some extended diffuse emission, while 
two galaxies of the same group (according to the spec- 
troscopy by Veilleux et al 2002) are within 90 kpc in pro- 
jection (Fig^J. Veilleux et al. (2002) classify this system 
as a post-merger, i.e. the tidal tails have been so diluted 
that they become hardly visible, while the center is still 
perturbed with a prominent knot of star formation. 

In this article, we adopt a standard flat cosmolog- 
ical model, with A = 0.7, and a Hubble constant of 
70kms -1 Mpc- 1 . IRAS 11582+3020 at z=0.223 is then 
at an angular distance of 743 Mpc, and V— 3.6 kpc. The 
luminosity distance is 1111 Mpc. 

3. Observations 

First, we performed IRAM 30m observations in May 2005 
(see Fig.|2). We used 2 SIS receivers to observe the 2 po- 
larizations simultaneously at 94.253 GHz, the redshifted 
frequency of the CO(1-0) line. At this frequency, the tele- 
scope half-power beam width is 27". The main-beam ef- 
ficiency is r? mb = TX/T mb =Q.77, and S/T* A = 6.1 Jy/K. 
The typical system temperature was 120 K (on the T£ 
scale). Wobbler switching mode was used, with reference 
positions offset by 2' in azimuth. Two 1 MHz filter banks 
provided a total bandwidth of 512 MHz, or 1600 kms -1 , 
with a velocity resolution of 3.2 kms -1 . 




Fig. 1. Red image of IRAS 11582+3020 from Kim et al 
(2002). The two objects marked "G" are galaxies in the 
same group (a fore-ground star is also indicated). The 
length of the horizontal bar at the bottom is 10"= 36 kpc 
at z=0.223. 
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Fig. 2. The IRAM 30m CO(1-0) spectrum (solid line), 
compared with the integrated spectrum from the interfer- 
ometer (dash line). The velocity is relative to z=0.223. 

IRAS 11582+3020 was subsequently observed with the 
6 antennae of the PdBI in January and February 2006 in 
the new A and B configurations. The 3mm receiver was 
tuned to 94.253 GHz and the 1mm to 240 GHz for the dust 
continuum (as no other CO line falls in the tuning bands). 
The dual-band SIS receivers yielded SSB receiver temper- 
atures around 40 K and 50 K at the two observed frequen- 
cies. The system temperatures were 120 K for CO(1-0) 
and 400 K at 240 GHz. Four correlator units covered a 
total bandwidth of 580 MHz at each frequency, provid- 
ing a nominal frequency resolution of 1.25 MHz (4kms _1 
for the CO (1-0) line), however we smoothed the signal 
to 28.6 kms -1 channels. On average, the residual atmo- 
spheric phase jitter was less than 20° at 3mm, consistent 
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Fig. 3. a) CO(l-O) integrated intensity contours in the inner ~7" (25.2 kpc) of IRAS 11582+3020. Contours are -0.5, 
0.5 to 4Jykms -1 beam -1 in steps of 0.5Jykms -1 beam -1 (l-er=0.2Jykms -1 beam -1 ). The black cross marks the 
position of the phase tracking center(aj2ooo=12' l 00 m 46.8 ;i and (5j2ooo=30°04'15.0"). The red cross marks the position 
of the CO peak intensity and CO-based dynamical center (a j20 oo = 12' l 00 m 46.85 s and 5j2ooo=30 04'14.3"). b) The 
mean velocity field derived from the 1-0 data using a 3<j-clipping in each channel map, in line contours, spanning the 
range (100 km s -1 , 500 km s -1 ) in steps of 100 km s -1 . The velocity scale is relative to z=0.223 with a CO-derived 
systemic velocity of v S j, s =300 km s -1 , i.e., zco=0-224 (redhifted by 180 kms -1 from the optical line which is at 
z=0.2234). The red line shows the position of the apparent major axis of the galaxy. Beam-sizes are represented by 
hatched ellipses. 



with a seeing disk of 0.34"-0.46" size and with a ~ 5% 
loss of efficiency. 

Reduction using the GILDAS software provided data 
cubes with 512x512 spatial pixels (0.14" /pixel) and 64 
velocity channels of 28.6 kms -1 width. The cubes were 
cleaned with the Clark (1980) method and restored by a 
1.3" x 1.0" Gaussian beam (with PA=54°) at 94.253 GHz 
and 0.43" x 0.25" (with PA=44°) at 240GHz. The mis 
noise levels in the cleaned maps (at 28.6 kms -1 velocity 
resolution) are 0.8 mJy beam -1 for the CO(1-0) line. No 
continuum emission was detected at 3mm (possible AGN) 
or at 1mm (possible dust emission) , down to rms noise lev- 
els of 0.1 mJybeam -1 and 0.5 mJybeam -1 in a 580 MHz 
bandwidth at 94 GHz and 240 GHz, respectively. Given 
the IRAS 100/im flux of 1.5 Jy, the upper limit at 1.2mm 
is compatible with a typical starburst SED. 

4. Results 

As revealed in FigEJ the total PdBI flux is ~ 50% of what 
is found with the single dish. When taking the intrinsic 
flux uncertainties into account this may indicate that some 
extended emission is not recovered by our interferometer 
observations. With a beam of 27" at the 30m, we measured 
a CO(1-0) flux toward the center of I(CO)= 1.3 K.kms -1 
(in the T^ scale). With a conversion factor of 6.1 Jy/K, 
the integrated flux is then S(CO)=7.8 Jy.kms -1 . With 
the PdBI, we measure 3.7 Jy.kms -1 in the FOV= 54" of 
the CO (1-0) map. This can be explained by the long PdBI 
baselines which are not sensitive to CO emission extending 



over 10" and detect only the more clumpy distribution. 
The optical structure does reveal the existence of such an 
extended component, and also hints at a weak diffuse tidal 
tail in this perturbed system, classified as a post-merger 
(Veilleux et al 2002). 

The CO luminosity for a high-z source is given by 

D 2 

L co = 23.5I C o^B n K km s -1 pc 2 

(1 + z) A 

where VLb is the area of the main beam in square arcsec- 
onds and Dl — \\\\ is the luminosity distance in Mpc. We 
compute H2 masses by Mh 2 = otLco M , with a = 0.9 for 
ULIRGs, instead of 4.6 for normal galaxies (i.e. Solomon 
et al. 1997). The molecular masses are listed in Tabled 

The integrated CO(1-0) map is plotted in Fig|3 to- 
gether with the isovelocity curves. The source is resolved 
at least in the direction of the beam minor axis, where the 
deconvoled size of the CO emitting region is of the order 
of 0.8" in diameter (~ 3 kpc). The signal-to- noise ratio 
is around 10 in most of the channels, implying a preci- 
sion on the position of the peaks of ~0.1" . Therefore it 
is possible to detect the shift of the barycenter in each 
channel map, with a kinematic major axis aligned at PA 
~135° (Fig|3|). The velocity gradient is also clearly seen 
in the position-velocity diagram taken along this position 
angle (Fig^J). Since ^135° is also the position angle of the 
extended optical isophotes (Fig [I}, all available data are 
compatible with a post-merger relaxed system. 

The CO integrated spectrum seems to show a double- 
horn profile, indicative of a rotating disk. Since the con- 
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Instrument 


AV 


Area 


M(H 2 ) 




kms -1 


Jy.km s _ 


M 


30m 


476±33 


7.8±0.5 


1.2±0.1 10 10 


PdBI 


550±52 


3.7±0.3 


6±0.4 10 9 



Table 1. CO results from IRAM-30m and PdBI, and deduced 
H 2 mass (AV is the FWHM). 
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Fig. 4. We show the CO(l-O) position-velocity diagram 
along the major axis of IRAS 11582+3020 (line contours 
range: 2.5a, to 10.5<r in steps of ler; 1<7=0.7 mJy beam -1 
in channels of 86 km s _1 ). Velocities have been re-scaled 
to zco=0-224 and angular offsets are relative to the dy- 
namical center. The white crosses mark the velocities cor- 
responding to the peak brightness temperatures as a func- 
tion of position. 

tours of the red image are suggesting an inclination of 
i=50° with the plane of the sky for a disk geometry, we 
can estimate a dynamical mass, assuming that the width 
of the CO profile is twice the projected rotational veloc- 
ity V rot sin i = 250 kms" 1 . With V rot = 326 kms -1 , 
inside a radius of 1.5 kpc, i.e. the extent of the central 
molecular disk, the indicative dynamical mass is Mdyn = 
3.4 10 10 Mq. With the low conversion ratio proposed for 
the ULIRGs, the gas mass in the central molecular disk 
(M(H 2 ) = 6 10 9 Mq) is a small fraction (~ 20%) of the 
dynamical mass inside 1.5 kpc. 

5. Discussion and conclusions 

The ~ 1" resolution CO observations of the ultra- 
luminous starburst IRAS 11582+3020 at z = 0.223 pro- 
vides interesting insight into its molecular gas content and 
morphology. The CO emission is resolved, spatially and 
kinematically: we see evidence for a rotating disk. In ad- 
dition to a central molecular disk of ^1.5 kpc in radius, 



there is tentative evidence for additional extended emis- 
sion over 10-15" as about 50% of the 30m flux are not 
recovered by the interferometer. Such an extended compo- 
nent could correspond to the perturbed optical morphol- 
ogy of the evolved merger. This source size contrasts with 
the 10 nearby ULIRGs, mapped by Downes & Solomon 
(1998) where the molecular gas is more concentrated, with 
typical radii of 300 to 800 pc. The derived total amount 
of molecular gas in IRAS 11582+3020 is comparable to 
the average mass in local ULIRGs of ~ 5 10 9 M Q and its 
Mgas/M^yra ratio is also comparable to the value of ~ 1/6 
found by Downes & Solomon (1998). 

The physical properties of IRAS 11582+3020 appear to 
be intermediate between that of local ULIRGs and high-z 
submillimeter galaxies (SMC) mapped by Tacconi et al 
(2006). The latter have also a typical radius of 2 kpc, and 
their FIR-to-CO luminosity ratio is higher than the local 
ULIRGs ratio, due to their higher FIR luminosities. There 
might be a trend, also followed by IRAS 11582+3020, 
for high-z galaxies to have a larger L(FIR)/M(H2) ratio 
(Riechers et al 2006), however this must be confirmed with 
larger samples. The high-z SMG of Tacconi et al (2006) are 
only mapped in the higher J rotation lines of CO, which 
may explain the complete absence of extended emission, 
even in the case of mergers. From their study of the three 
CO-brightest z>4 QSOs, Riechers et al (2006) limit a po- 
tential extended CO emission to <30% of the total. Two 
of their quasars are gravitationally lensed (with magnifi- 
cation factors of ~3-7), i.e. differential magnification may 
slightly increase this limit. 
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